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Abstract 

We develop the cluster self-consistent field method incorporating both electronic and lattice degrees of 
freedom to study the origin of ferromagnetism in Cs2AgF4. After self-consistently determining the har- 
monic and anharmonic Jahn-Teller distortions, we show that the anharmonic distortion stabilizes the stag- 
gered x 2 -z 2 /y 2 -z 2 orbital and ferromagnetic ground state, rather than the antiferromagnetic one. The am- 
plitudes of lattice distortions, Q2 and Q3, the magnetic coupling strengthes, J A V , and the magnetic moment, 
are in good agreement with the experimental observation. 

PACS numbers: 71.70.Ej,75.10.-b,75.30.Et 



1 



I. INTRODUCTION 



Recently, layered perovskite compound Cs2AgF 4 containing the spin- 1/2 4d 9 Ag(II) ion has 
attracted great interest, owing to its isostructural to the high-T c 3d 9 Cu(II) cuprates. Cs 2 AgF 4 was 
first refined to the tetragonal structure with the space group I4/mmm in 1974 by Odenthal et al. yj]. 
Very recently, however, McLain et al. [2] found that the crystal structure is orthorhombic with the 
space group Bbcm. The magnetic susceptibility and the inelastic neutron scattering experiments 
showed that this compound is well described as a two-dimensional (2D) ferromagnet (FM) with T c 
~ 15 K yfl. Further, McLain et al. suggested that Cs 2 AgF 4 is orbitally ordered at all temperatures 
of measurement, and the orbital ordering (00) is responsible for the FM. These properties is in 
sharp contrast to the high-T c parent La 2 Cu0 4 that are antiferromagnetic (AFM) insulators. The 
microscopic origin of the unusual FM and 00 in Cs 2 AgF 4 attracts a lot of attentions. 

In Cs 2 AgF 4 , the basal plane consists of a 2D lattice of Jahn-Teller (JT) distorted AgF 6 octahedra 
with a pattern of alternating short and long Ag-F bonds. These analogous compounds with orbital 
degeneracy, such as Cs 2 AgF 4 and K 2 CuF 4 , etc. turn out to be 2D FM [2], while other compounds, 
such as K 2 NiF 4 and Rb 2 MnF 4 , etc. are AFM Q3|, [4£| . Notice that the former with active JT ions leads 
to an orthorhombic structure, while the latter with nondegenerate orbitals only has a tetragonal 
structure. Although the difference between these two structures seems small, i.e. the main change 
is the position of the fluorine atoms in the basal plane, such difference is of the utmost importance 
in leading to the distinct properties |5J. From early studies [(1 7, 8, 9] in perovskite compounds, it 
is known that the anharmonic JT effect is a decisive factor for the orthorhombic crystal structure. 
To date, a lot of studies have been done for Cs 2 AgF 4 utilizing the density functional theory lioi 
[hIQQQ[l5|l however, few authors focus on the anharmonic JT effect on the lattice, 00 and 
FM groundstate properties in Cs 2 AgF 4 . 

Previous studies in analogous compound K 2 CuF 4 have provided us two distinct scenarios for 
the FM and 00 ground state. On the one hand, Kugel and Khomskii showed that the cooperative 
JT effect, especially the anharmonic JT effect Jlfj], play a key role in stabilizing the 00 ground 
state in K 2 CuF 4 . On the other hand, to resolve the difficulty of the Kugel-Khomskii (KK) elec- 
tronic superexchange (SE) coupling model [17Q, which is usually suitable for the Mott-Hubbard 
insulators with \J d < A (here is the on-site Coulomb repulsion and A is the charge transfer 
energy), and to address the FM and 00 ground state, Mostovoy and Khomskii fl 1 311 proposed a 
modified SE coupling model, which is suitable for the charge transfer insulating K 2 CuF 4 with 
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> A. Obviously, to understand the unusual groundstate properties, such as the lattice, orbital and 



magnetic properties in K 2 CuF 4 and Cs 2 AgF 4 , one should incorporate the cooperative JT effect Ill9h 
and the anharmonic JT effect into the charge-transfer SE interactions. 

To treat these strongly correlated systems more precisely, based on the cluster self-consistent 
field approach developed previously [20], we explicitly take into account the orthorhombic JT dis- 
tortions and the charge-transfer SE interactions, in which spin order, 00 and lattice distortion are 
determined self-consistently. We show that driven by strong anharmonic effect and Hund's cou- 
pling, Cs 2 AgF 4 has a much more stable ferromagnetic ground state. The theoretical amplitudes of 
lattice distortions, Q 2 and Q 3 , the magnetic coupling strengthes, J x>y , and the magnetic moment, 
are in good agreement with the experimental observation in Cs 2 AgF 4 . This paper is organized 
as follows: an effective Hamiltonian and the cluster self-consistent field (SCF) approach are de- 
scribed in Sec. II; then the results and discussions about the lattice structure and distortions, orbital 
ordering and magnetic properties in Cs 2 AgF 4 are presented in Sec. Ill; the last section is devoted 
to the remarks and summary. 

H. MODEL HAMILTONIAN AND METHOD 

An effective 2D Hamiltonian in Cs 2 AgF 4 including both the SE and the JT couplings is written: 

H = Hqj T + H SE (1) 

where H SE describes the highly symmetrical effective SE couplings between spins and orbitals; 
Hqjt describes the JT couplings associated with the electron and the lattice distortion. Note that 
spin, orbital and lattice degrees of freedom couple to each other in Hamiltonian (1). The JT effect 



associated with one hole per Ag i+ site in the AgF 4 sheets reads 1190 



Hqjt = g £(fiaT? + Q a <i*) + | + Q%) 

i i 

+G £[(g| - Ql)^ - IQaQirf] (2) 



where both the linear and the quadratic vibronic coupling terms have been included. The first and 
second terms describe the linear harmonic JT effect. The third one, i.e., the quadratic coupling, 
arises from the anharmonic JT effect and contributes to the anisotropic energy [1(3, [ijj]. Here 
Q, 2 and Q, 3 are the normal vibration coordinates, defined as Q 2 = {—X\ + X 2 + Y 3 - F 4 )/2 and 
g 3 = (-Xj +X 2 -Y 3 + Y 4 + 2Z 5 - 2Z 6 )/ Vl2 UM with X, Y and Z being the coordinates of the 



i - th F ions. And g is the linear JT coupling strength, G is the coefficient of quadratic coupling, 
and K is the elastic constant. For Cs2AgF 4 , we fix K = 10 eV/A 2 throughout this paper. 

Since Cs 2 AgF 4 is strongly correlated charge-transfer insulator [11], the 2D SE interactions 



between 4d orbitals of Ag 2+ ions through 2p orbitals of F atoms are described Ill8ll : 

Hse = 2 [(Ji + hit + hl-ll a )$i ■ $ 



' r+a 



i,a 
a=x,y 



+/ 4 /f + / 5 /f/)U (3) 

here the operator 5*, denotes the S = 1/2 spin at site i, and If = cos (2nm a /3) -f. - sin (2nm a /3) t*, 
which is the combination of the components of the orbital operators t. a = x,y, and (m x , m v )=(l, 
2), denote the direction of a bond in the AgF 4 sheets. In Eq. (3), the coupling coefficients 
J„ (n= 1-5), read: J, = t 2 [l/U d + 2/(2A +U P )- J H /(2U 2 d )], J 2 = 4t 2 [l/U d + 2/(2A + U p )\ 
h = 4t 2 [l/U d + 2/(2A+ U p ) + J H /(2U 2 )}, J 4 = f 2 [2/A-l/f/ J -2/(2A+C/ p )] and 7 5 = 
t 2 [l/U d + 2/A - 2/(2A + U p ) + 3J H /(2Uj)] with t=t 2 pd /A, respectively. Here A is the electron 
transfer energy between silver 4d to fluorin 2p orbitals; U d , \J P are the Coulomb repulsion ener- 
gies on Ag and F, respectively; and J# is the Hund's coupling on Ag. Note that the coefficients 
may be positive and negative, thus the AFM and the FM components coexist in J„. Since the 4d 
orbitals of Ag 2+ ions are less localized than the 3d orbitals of Cu 2+ ions in K 2 CuF 4 , the Coulomb 
repulsion LL of Ag 2+ ion is about 3 ~ 5 eV, and Jh ranges from 0.1 eV to 0.5 eV, in accordance with 

nrr 



the LDA+U results llOLlll, 



15J]. In general, the Coulomb repulsion XJ P is about 5 eV in F . The 



charge transfer energy A is roughly estimated from the difference of the centers of gravity between 
the 4d and 2p levels in band structures, which also allows to estimate t pd from the bandwidth of 
electronic structures available 1150, A=U rf -2 eV, and t pd - 0.6 eV. 

In such a strongly correlated spin-orbital-lattice system, according to Feynman-Hellman theo- 
rem, the ground state energy is minimized with respect to Q,2 and Q, 3 , i.e., 

<^ >= 0.<|f>=0. 
dQn oQ B 

From which one could find the strength of the normal modes and the lattice distortion critically 
depend on the orbital polarization through the following equations 

KiT^ + AGir^ir]) 



(Qa) 



' 4G 2 ((t*) 2 + <rp 2 ) - K 2 



K{t]) + 2G((r-) 2 - (r^) 2 ) 

g 4G 2 (<T * >2 + (T z )2) _ K 2 ^ 



Notice that in the absence of the anharmonic JT effect (G=0), <Q,2> ~ -(g/K)<r ; v > and <Q,3> ~ 
-(g/K)<n>. To obtain the amplitude of the lattice distortion, one should determine the spin, the 
orbital and the deformation configurations self-consistently. 

In order to treat the spin-orbit correlations and fluctuations with high accuracy in strongly 



correlated system, the cluster-SCF approach Il20ll developed previously is applied to the spin 



orbital Hamiltonian (1). The cluster-SCF approach includes the exact treatment to the interactions 
inside the cluster, and the self-consistent field treatment on the interactions between the cluster and 
the surrounding environment. The main procedure is shortly outlined as follows: first, we choose 
a cluster consisting of 4 Ag 2+ ions. The i - th site is surrounded by two Ag sites inside the cluster 
and two Ag sites outside the cluster, which provide two internal interactions and two external SCF 
fields, respectively. Thus, the effective Hamiltonian of the cluster reads: 

Kluster = ^[(Jl + jV? + + Jtflf+atfi-ZiH, 

i,a ~ 

a=x,y 

+J 4 If + J 5 IfIf +a ] + g Yj.Qirt + Qht*) 

i 

+G^KQ? 3 -QlX-2Q i2 QBT*l 

i 

+§ + g 2 *) + 2>r' (5) 

i i 

with the SCF contributing from the / - th external site interacting with the i - th internal site 
through Hi/, h. = Try{pyHij>), where p/ denotes the reduced density matrix of the j'th site, and 
i runs over all sites inside the cluster. We first substitute the spin coupling S; • Sj +a into the cluster 
Hamiltonian with the initial spin correlation function (Si • S*i +a ), and diagonalize the orbital part of 
the cluster Hamiltonian (5) in the presence of the orbital SCF. The orbitalization (f) and the orbital 
correlation functions (r, • r i+a ) are thus obtained. And then, we substitute the orbital operator and 
the orbital couplings with (r) and (f, • f i+a ) into h duster , and diagonalize the spin part of the cluster 
Hamiltonian in the presence of the spin SCF. Hence, we obtain a set of new averaged spin 5, and 
spin correlation functions ^r^i +a . Repeat the above steps until the groundstate energy and the spin 
and orbital correlation functions converge to the accuracies. From the stable spin-orbital ground 
state, one could get the spin coupling strengths J* and J y through J a = Ji+J 2 < If >+ h< Iflf +a >> 
here a = x, y. The advantage of our approach over the traditional mean-field method is that the 
short-range spin and orbital correlations as well as the quantum fluctuations are taken into account 
properly, especially in the low-dimensional systems. 



III. RESULTS AND DISCUSSIONS 



In this section, we present the numerical results of the ground state of the charge-transfer in- 
sulator Cs 2 AgF 4 within the cluster-SCF approach. We mainly discuss the role of the anharmonic 
effect on lattice structure and distortions, orbital ordering and magnetic properties in the ground 
state. 



A. Groundstate Phase Diagram 

To clarify the roles of the SE coupling and the JT effect on the origin of the FM, we first perform 
numerical calculations on Eq. (5) and obtain the g-J# phase diagram, as shown in Fig. 1 . Since 
the on-site Coulomb interaction in Cs 2 AgF 4 is not well defined, we present the g-J# phase 
diagram for U</=3, 4, and 5 eV in Fig. 1(a), 1(b) and 1(c), respectively. Accordingly, we find that 
the anharmonic distortion strength G plays a very important role in the phase diagram. With the 
variation of G from negative to positive via zero, the FM-AFM phase boundaries exhibit critical 
changes. As we see in Fig. 1(a), at G=-0.75g, the stable groundstate phase in the regions I, II, III 
and IV is AFM, while that in the region V is FM; at G=0, a vertical line separates the AFM ground 
state in the regions I, II, and III from the FM one in the regions IV and V; when the anharmonic 
distortion strength G becomes positive, more regions become FM ordering; at G=0.75g, only the 
ground state in the region I is AFM. Fig. 1(b) and Fig. 1(c) qualitatively resemble to Fig. 1(a). 

The competitive FM and AFM couplings in the SE Hamiltonian result in the FM phase for large 
J# and the AFM phase for small J H in these phases. As we expect in Fig. 1(a) to Fig. 1(c), strong 
Coulomb repulsion favors AFM phase, hence the AFM regime becomes large with the increase of 
the U d . Besides the fact that strong Hund's coupling J H favors the FM phase, the large anharmonic 
JT effect also favors the FM ground state. Therefore, the complete consideration of the realistic 
interactions in Eq. (1) most probably leads to the FM ground state in Cs 2 AgF 4 . And from the 
recent first-principles electronic calculations, the Coulomb repulsion between Ag 4d electrons is 



about 3 eV fl lOfl . We expect that the Hund's coupling J H lies in 0.1 ~ 0.3 eV, so that the interaction 
parameters of Cs 2 AgF 4 fall in the FM region in the phase diagram. 

Experimentally, it is hard to determine the sign of the anharmonic coupling strength G in 
Cs 2 AgF 4 . From Fig. 1, it is seen that when positive anharmonic effect is very large, the ground 
state stabilizes in the FM phase even in the absence of Hund's coupling. Therefore the positive 
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FIG. 1: (Color online) Phase diagrams of JT coupling g versus Hund's coupling J H in Cs 2 AgF4 for different 
Coulomb repulsions (a) U</=3 eV, (b) Urf=4 eV, and (c) Uj=5 eV. The solid circles, dashed line, solid 
triangles and solid squares denote the FM/AFM phase boundaries for the anharmonic coupling G=-0.75g, 
0, 0.25g and 0.75g, respectively. G is in the unit of eV/A 2 . 



anharmonic coupling most favors the FM ground state. On the contrary, for G < 0, due to the SE 
coupling, the ground state is FM when the harmonic JT coupling strength g is small; however, 
when g becomes large, the anharmonic coupling favors the AFM ground state. Thus, considering 
the strong JT effect in Cs 2 AgF 4 , one expects that G is most probably positive, which will be fur- 
ther confirmed in the lattice distortion, the orbital ordering and the magnetic properties in what 
follows. 



B. Lattice Structure and Distortions 

Different from La 2 Cu0 4 , the unusual FM and the orthorhombic structure of Cs 2 AgF 4 is the 
result of the interplays of the spin, orbital and lattice degrees of freedom, especially the contri- 
bution of the anharmonic JT effect. The normal coordinates Q 2 and Q 3 of the JT distortions in 
Cs 2 AgF 4 are obtained self-consistently through Eq. (4) and (5). We find that for positive and large 
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G, the structural distortion in the ground state is antiferro-type, (Q 2 , Q3; -Q2, Q3), corresponding 
to the alternative short and long bonds in the AgF 4 plaquette. However, the groundstate energy 
of the ferro-type (Q 2 , Q3; Q2, Q3) distortion is higher than that of the antiferro-type distortion, as 
shown in Fig. 2. And the ferro-type distortion only appears in the strong JT region because of its 
instability in this system. 




g(eWA) 



FIG. 2: (Color online) Dependence of groundstate energy on the JT coupling strength with G = 0.75g in the 
antiferro-type (a) and the ferro-type (b) distortions in the AgF 4 plaquette. The dashed vertical line indicates 
the appearance of the ferro-type distortion. Theoretical parameters: U=3 eV and J#=0.3 eV. 



Notice that the experimentally observed lattice structure has the antiferro-type distortion with 
the long and short Ag-F bond lengths of 2.441 and 2.111 A in the AgF 4 plaquette [2], one can 
deduce the corresponding distortions Q 2 « 0.33 A and Q 3 « -0.19 A. Theoretically, through the 
self-consistent numerical calculations, we obtain the amplitudes of the Q 2 and the Q 3 distortions 
in Cs 2 AgF 4 with the variation of the JT coupling for different anharmonic coupling parameters, 
as shown in Fig. 3. When the linear JT coupling strength g is about 4.8 eV/A at G=0.75g, the 
theoretical distortions are |Q 2 | ~ 0.33 A and Q 3 « -0.19 A, see the dashed vertical line in Fig. 3, 
in good agreement with the corresponding experimental observation in Cs 2 AgF 4 [2]. Meanwhile, 
we find that the positive G leads to the negative Q 3 , giving rise to the correct distortion c/a = 
2s/(l+s) < 1 in orthorhombic Cs 2 AgF 4 . However, the negative G gives a positive Q 3 , resulting 
a wrong distortion c/a = 21/(l+s) > 1. All the other analogous substances without JT ions, e.g. 
K 2 NiF 4 , K 2 MnF 4 Q3L |4|], are tetragonal, rather than the orthorhombic symmetry in Cs 2 AgF 4 with 
JT ions. In fact, it is the positive anharmonic effect that lowers the lattice symmetry, leading to the 
compression of the ligand octahedron along the c-axis [|22|l . 
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FIG. 3: (Color online) The amplitudes of IQ2I and Q3 as a function of the JT coupling with different G. 
Two dashed horizontal lines correspond to the experimental values of IQ2I and Q3; the vertical dashed line 
indicates the linear JT coupling g=4.8 eV/A for G=0.75g in Cs 2 AgF 4 . Other parameters are the same to 
Fig. 2 

C. Orbital Ordering 

The unusual FM and distinct orthorhombic structure in Cs 2 AgF 4 are in fact associated with the 



formation of the 
calculations 



ong-range 00, as shown by a few authors recently utilizing the first-principles 



tne long-rai 



1511 . Our study also confirms this point. More further, we demonstrate 
that the anharmonic effect plays an essential role in the 00 properties of Cs2AgF 4 . In general, one 
can describe the orbital occupied state of each site in terms of |0)=cosf |3z 2 - r 2 ) + sin||x 2 - y 2 ) 
with the orbital angle 0. Here '+' refers to the two sublattices of the antiferro-distortion in AgF 4 
plaquette. Fig. 4 shows the orbital angle 6 as a function of the linear JT coupling strength g with 
different anharmonic coupling G. The role of the anharmonic coupling on the orbital angle, hence 
the 00, is clearly seen. For strong JT distortion and G=0.75 g, the 00 is staggered z 2 - x 2 /z 2 -y 2 , 
see inset (b) in Fig. 4. The orbital angle obtained is about 61.8° for g=4.8 eV/A, as indicated by the 
dashed vertical line in Fig. 4, which is consistent with the experimental suggestion for Cs 2 AgF 4 
|2D. This is different from the case in La 2 Cu0 4 , where the distance between apical O and Cu is 
much larger than that in the a-b plane, resulting in the pure 3z 2 - r 2 hole orbitals. 

As a comparison, for G < 0, the orbital angle 9 > n/2; as |G| becomes very large, 6 — > 2tt/3, 
and the holes orderedly occupy the d 3A -2_ r 2/d 3v 2_ r 2 orbitals, as seen the inset (a) in Fig. 4. For G = 0, 
however, the orbital angle is nearly constant with the increase of the linear JT coupling g, giving 
rise to 0=68°. Note that in the absence of the anharmonic effect, the SE couplings contribute to the 
same 00 as the harmonic JT effect does. In this situation, the orbital configuration is characterized 
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FIG. 4: Orbital angle as a function of the JT coupling with different anharmonic coupling. Inset (a) and 
(b) correspond to the 3x 2 - r 2 /3v 2 - r 2 and x 2 - z 2 /y 2 - z 2 orbital orders for large negative and positive G 
values, respectively. The vertical dashed line corresponds to the case of g=4.8 eV/A for Cs2AgF4. Other 
parameters are the same to Fig. 2. 

by orbital angle 6 : tan6 = Qq/Qb [19]. Then QalQa ~ T J/ T r wn i cn is almost independent of 
the linear coupling g, indicating that the linear JT coupling has a little influence on the orbital 
angle in the absence of the anharmonic effect. This also shows the key role of the anharmonic 
effect in the 00 from the other aspect. 



D. Magnetic Properties 

The magnetic properties are also strongly influenced by the anharmonic effect, as shown in Fig. 
5. In the present situation, we find that the difference between J x and J y is negligible. From Fig. 5, 
one finds that in the absence of the anharmonic effect (G=0), due to the SE coupling, the magnetic 
couplings between Ag spins are FM, and nearly remain unchange with the increasing of g. While 
for G < 0, due to the SE coupling, the spin coupling between Ag ions is FM at small g. With the 
increase of the JT effect, the magnetic couplings between Ag ions, J x<y , vary from FM to AFM, 
which is in contradiction with Cs 2 AgF 4 . For g=4.8 eV/A and G=0.75g, as denoted by the dashed 
vertical line in Fig. 5, the magnetic couplings in the AgF 4 plaquette, i.e. J x v , are about -3.3 meV, 
very close to the experimental values, -3.793 — 5.0 meV Q2Q . The agreement between the present 
spin coupling strengths and the experimental observation is much better than those inferred from 



the first-principles electronic structure calculations 113 



This result, combining those of the lattice distortion and the orbital angle , is in accordance with 
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FIG. 5: (Color online) Magnetic coupling strengthes J T} . in the AgF4 plaquette as a function of the JT 
coupling with different anharmonic coupling. The vertical dashed line corresponds to the case of g=4.8 
eV/A for Cs2AgF4. Other parameters are the same to Fig. 2. 



the Goodenough-Kanamori- Anderson rules [24]. Also in Fig. 5, one notices that in the presence 
of strong anharmonic coupling G, the FM coupling strength increases with the linear JT coupling 
g. Therefore the FM couplings in Cs 2 AgF 4 further verify the key role of the anharmonic effect 
in the groundstate properties. In addition, the calculated magnetic moment is about 1 fi B , which 
approaches the classical value and is comparable to the experimental data observed in muon-spin 



relaxation experiment Ii25|l . This is a classical value of Ag ion, in comparison with the experimen- 
tal value 0.8 jx B at 5 K [2D. The reduction of the experimental magnetic moment may originate 
from three reasons: the presence of non-magnetic impurity phases and the co valence effects as 
pointed by Mclain, et al. in Ref. [2], as well as the AFM fluctuations between AgF 4 layers. We 
expect that the weak interlayer AFM interaction and the covalency effect, which are neglected in 
the present 2D model, contribute to the major part of the reduction of the magnetic moment of Ag 
spins. 



IV. REMARKS AND SUMMARY 



As we have shown in the preceding sections, the anharmonic JT effect drives the 4d orbitals 
of Cs2AgF 4 into the d z 2_ x i/d z 2_ y 2 ordering, rather than the d 3z 2_ £ 2 orbital ordering in La 2 Cu0 4 . In 
consistent with the Goodenough-Kanamori- Anderson rules [24], the ground state of Cs2AgF 4 is 
unusual FM insulating with the orthorhombic structure, different from the Neel AFM insulating 
ground state with the tetrahedral structure in La 2 Cu0 4 . Our results preclude the possibility of the 



11 



FM ground state originating from the co valency effect Ill2h . One consequence of such a difference 
is that the electron or hole doping in La 2 Cu04 leads to the strong AFM fluctuations, which may 
contribute to the Cooper-pairing glue for the high-T c superconductivity. We anticipate that the 
electron doping in Cs 2 AgF 4 will lead to weak FM fluctuation, though it will not contribute to the 
same SC mechanism as the cuprate superconductors. 

Notice that in the present intermediate correlated and charge transfer insulator, the 2D effec- 
tive SE interactions incorporating the JT couplings underestimate the covalency effect between 
Ag 4d and F 2p orbitals. The covalency effect between Agand F ions can be well considered 



15I] . On the other hand, 



within the first-principles electronic structure calculations |10L I12L 
once the JT distortion occurs, the hopping integrals, t x iV , between Ag ions are slightly different in 
different crystallographic axes. We expect that this change will not significantly modify our results 
quantitatively. 

In summary, an effective model Hamiltonian with spin, orbital and lattice degrees of freedom 
coupling to each other, allows us to self-consistently describe the lattice structure change, orbital 
ordering and magnetic coupling properties, etc. Driven by the JT distortion, especially by the 
anharmonic effect, Cs2AgF 4 stabilizes in the 2D FM ground state. The correct lattice structure, or- 
bital ordering and magnetic coupling observed experimentally can also be addressed in the present 
theoretical framework. 
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